Excitation spectroscopy of vortex lattices in a rotating Bose-Einstein condensate 
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Excitation spectroscopy of vortex lattices in rotating Bose-Einstein condensates is described. We 
numerically obtain the Bogoliubov-deGenne quasiparticle excitations for a broad range of energies 
and analyze them in the context of the complex dynamics of the system. Our work is carried out 
in a regime in which standard hydrodynamic assumptions do not hold, and includes features not 
readily contained within existing treatments. 
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The collective excitations of many particle systems 
play a central role in the understanding of the system's 
bulk properties and provide a basis for interpreting com- 
plex microscopic dynamics. In the case of a "super" - 
system - a superfluid or a superconductor - studies of 
excitations can be used to probe particle-particle inter- 
actions and afford essential insight into energy flow and 
dissipation. Furthermore, the novel topological character 
of excitations such as vorticies are of intrinsic fundamen- 
tal interest. Stimulated by the achievement of vortex lat- 
tices in trapped atomic Bose-Einstein condensate (BEC) 



I, there is an intense effort to understand quantized 
vortices and their excitations in this system |3| . Al- 
though vortex lattices have been extensively studied in 
superfluid ^He and type-II superconductors, most treat- 
ments rely on approximations that do not adequately re- 
flect the situation of the trapped atomic BEC. For exam- 
ple, in the atomic BEC, the motion of the vorticies (and 
vortex lines) cannot be completely isolated from the bulk 
motion of the fluid in which they are imbedded. To ad- 
dress these limitations, a number of studies have been 
presented Jp, Q, including treatments of Tkachenko os- 
cillations [g and hydrodynamic shape oscillations [9lllo|. 
Nevertheless, a full analytical solution, valid over a wide 
range of energies and densities, is difficult if not impos- 
sible to construct due to the inherent complexity of the 
system. In this paper we describe a comprehensive ap- 
proach based on a numerical calculation. We concentrate 
on the regime in which various hydrodynamic assump- 
tions, such as the rigid body approximation, break down. 
Our method enables us to search for new effects and to 
consider the case of a small condensate with few vortices 
of varying core size. 

We first present a general description of the normal 
modes based on the observed structure of the quasipar- 
ticles. We then discuss the excitation spectrum. We find 
that the absence of the rigid body rotation of the fluid 
modifies the symmetry of the energy spectrum. We also 
find that the spectrum is decorated by features arising 
from resonant processes that occur when wavelength of 
the surface excitations matches the intervortex distance. 

Since we are interested in the rotating condensate, a 
rotating reference frame is introduced and so the grand 



canonical free Hamiltonian F includes additional cen- 
trifugal term: 

F = J dr^;^n-^l- no + ^gi!^^ip, (I) 

where H = -{h^/2M)\7^ + Vtrir),C, = -ihd/dcl), n 
is the rotation frequency of the reference frame, g — 
inah^ /M is the interparticle coupling constant with the 
s-wave scattering length a, ip is a field operator satisfy- 
ing the bosonic commutation relation and Vtr(r) is the 
trapping potential. We use harmonic trapping potential 
Vtr = MujIX^"^ + ?/^)/2 + Mu;'^z'^/2. Using the standard 
procedure in which the mean field ■0(r) = ( '0(r, t) ) and 
the fluctuation (j){r,t) = ip{r,t) — ipir) are defined, the 
Gross-Pitaevskii equation (CPE) in the rotating frame 
at T = is given from the free energy [Eq. |^] to be 



(2) 



while applying the Bogoliubov transformation, (f> (r, t) = 
X] Ui(i")aje~*'^j* — Uj*(r)aj^e*"j*, yields coupled Bogoli- 
ubov equations in the rotating frame. 



{H + 2gijj*ifi — /i — nCz) Uj — gi'^Vj 
{H + 2gifi*ip — /i + nCz) Vj — gip*'^Uj 



hujjUj 
—hujjVj. 



(3) 



These equations, written as sparse matrices of order 
10^ X 10^, are solved numerically using parallel computing 
packages |0. All calculations are done in a quasi 2- 
dimensions (2D) limit consistent with the experimentally 
relevant pancake geometry All values given in this 
paper are based on the dimensions of the trap potential: 
y^h/MuJti for the length and w^"^ for time. For these 
2D calculations, the effective coupling constant (72D is 
used, which is defined 52D = g{mujz/2'KTiY/'^ 0. We 
use Wtr = 27r X 74.25IIz, oj^ = lOwtr and 10^ rubidium 
atoms with a — 5.77nm which gives 520 = 0.0724 in trap 
potential dimensions. 

To calculate the ground state we propagate an initial 
trial state in imaginary time, a technique which con- 
sistently decreases the energy until convergence is ob- 
served. For the above parameters, a triangular vortex 



lattice state containing 31 vortices with a six- fold sym- 
metry is found for a trap rotation frequency fl = 0.8, 
while one with 7 vortices with the same symmetry is 
found for = 0.5. The lattice with 31 vortices has a 
Thomas- Fermi (TF) radius ttf = 7.1, a healing length 
^ = 0.23 with an intervortex spacing 6 = 2.1; correspond- 
ing values for the lattice with 7 vortices are ttf =5.9, 
^ = 0.19 and b = 2.7. It is clear from the relatively large 
ratio ^/b that these systems are not strictly in a hydro- 
dynamic regime and, especially for the case of 7 vortices, 
that the large b/rxF ratio indicates that the system does 
not have sufhcient number of vortices to ensure a rigid 
body rotation of the fluid. 

We can visualize the evolution of the density profile of 
the jth normal mode using the time-dependent probabil- 
ity density IV' + w^e"*'^'* — u*e''^^*|^, where V'(r) indicates 
the ground state wave function and ujj represents the 
corresponding eigenenergy. Examining the time evolu- 
tion, we observe that the Tkachenko modes occur only 
in the lower energy regime. The time dependent density 
profiles for the Tkachenko modes were found to be in 
excellent agreement with the experimental results from 
JILA 01 and are described elsewhere 0. We also find 
that the hydrodynamic shape excitations which involve 
significant motion of the bulk fiuid are generally found 
throughout the entire energy range. 

Following we note that in the absence of the lat- 
tice, hydrodynamic bulk excitations may be classified by 
the radial and angular quantum numbers n and m. In 
the presence of the vortex lattice we seek a similar de- 
scription. First, we find that the amplitudes for u and 
V, when plotted as a function of position, typically show 
an annular symmetry, with density dips at the position 
of the individual vortices. In the radial direction, the 
nodes in u and in v are quite clear and the number of 
nodes can be used to define the radial quantum numbers 
n. We next observe that the phase of u is no longer the 
same as that of v. Specifically, we find that if we write 
tp = ipoe^^'" then u and v are uoe*'^+''^° and voe'^^~^^''' re- 
spectively, where -ipo, U(uVo, S and Sq are real functions 
of r; n' = ij}*u — ipv |l| is therefore adopted to char- 
acterize the phase properties of the excitations. Phase 
plots of n' generally indicate that a well-defined multi- 
ple of 27r-phase change occurs around the edge of the 
cloud, even though the angular dependence of n' is often 
more complicated than e""*^. Mathematically, m may be 
defined by a contour integration of V5 along a closed 
contour enclosing the condensate at the Thomas-Fermi 
radius: m = ^ §f.dT ■ WS. Temporal evolution of hydro- 
dynamic shape modes have shown that they correspond 
well to circulating (annular) waves with \m\ wavelengths 
and that the sign of m determines the direction of the 
rotating wave; they rotate counter-clockwise if m > 0, 
and clockwise if to < 0. 

The simplest examples of the hydrodynamic excita- 
tions are the dipole (or the "center-of-mass" ) modes 
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FIG. 1: (a) Dipole mode (n = 0, m = 1) with lo = l-fi = 0.20 
and (b) quadrupole mode (n = 0, m = 2) with cu = 0.36. 
Surface plot of phase from to 2n (dark to light) across n' 
(c) for the dipole mode and (d) for the quadrupole mode. It 
is clear that the phase change along the circle drawn in each 
plot is 2n for m = 1 and 47r for m — 2. 

(n = 0, m — ±1) and the quadrupole modes {n = 0, 
TO = ±2) [13]. The TO = -1-1 dipole mode displays two 
peaks in condensate density which rotate at a frequency 
uj = 1 — Q while the vortices remain fixed. Meanwhile 
for the TO = — 1 mode, the whole condensate "sloshes" 
around a small circle together with the vortex lattice 
(without any visible lattice deformation) and occurs at a 
frequency lu = 1 + Q. In the lab frame, the frequencies 
for the two dipole modes are simply the trap frequency, 
which is 1 in our dimensions. Figures ^^a) and (b) show 
an instantaneous density profile ( "snapshot" ) taken from 
the temporal evolution of the to = -1-1 dipole and the 
TO = +2 quadrupole modes respectively. Figures ^c) 
and (d) depict the corresponding phases of the noncon- 
densate fiuctuation density n' . Along a circular contour 
drawn in each plot at Thomas- Fermi radius, an integral 
number of 27r-phase-changes can be equated to to. The 
{n = 1, TO = 0) and (n = 2, to = 0) breathing modes were 
also be identified and occurred at to — 2 and lu — 2.7 re- 
spectively (for fl = 0.8). For the breathing mode, the 
vortices move sympathetically with the bulk fluid mo- 
tion, and there is no difference between the frequencies 
in the rotating and nonrotating reference frames. 

We note additionally that a surface mode with quan- 
tum number ±to introduces \m\ nodes on the edge of the 
condensate, each of which involves a phase singularity 
with a ±27r-phase change around that singularity. Al- 
though these nodes cannot be identified as vortices, par- 
ticularly for small to, we find that as to increase, we see 
TO vortices gradually forming. We conjecture that this 
effect is linked to the vortex nucleation and annhilation 
processes, particularly when a change of rotational speed 
is imposed on the system. For sufficiently large to, the 
vortices merge with each other in a continuous manner 
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FIG. 2: Density plot for n = and (a) m = 30 (b) m = 50 
mode, (a) Surface wave forms visible vortices, (b) Vortices 
around condensate merge into each other. 



and appear to form a vortex sheet (Fig. |2Jl 14]. How- 
ever, no clear evidence of the phase transition to vortex 
sheets has been identified. 

We construct the excitation spectrum in the form of 
a plot of the excitation mode angular momentum as 
a function of energy - the central result of this pa- 
per. The angular momentum of the excitations with 
respect to the rotating condensate may be defined 
as = [{lu - l^)Nu + {h + k,)N^] I {Nu + N.^), where 
la — (/ dr a*{—ihd/d(j))aj j (J dv a*a) with a = u,v,ip 
and Na — Jdrjup, A'^u = Jdrjwp. This gives the average 
difference in the angular momenta for the quasiparticles 
and the condensate ground state. captures the effect 
of the detailed vortex lattice structure while m is an effec- 
tively coarse-grain-averaged result. The angular momen- 
tum Lz vs. excitation energy Wn.m curves for Q = 0.5 (7 
vortices) and H. — 0.8 (31 vortices) and the corresponding 
m vs. LUn^m curves in the rotating frame are presented 
in Fig. 121 The points with the same radial quantum 
number n are connected such that different points on 
each line correspond to the same n but to different to. 
The leftmost curve corresponds to the surface mode with 
n — 0, while the other curves are modes with higher n. 
In addition to these curves, we identify the Tkachenko 
modes as the group of points (the open squares) in the 
bounded region of very low energy. The points which are 
not connected (open triangles) are those modes which 
display very strong coupling of the bulk fluid motion to 
the vortices, thus making exact identification of the mode 
character difficult. This happens particularly for small m 
with higher n modes, for which excitation is concentrated 
deep within the condensate. 

In our previous work we found that, in the rotat- 
ing frame, there is a small energy splitting 20 between 
the +m and —m states independent of |to|, resulting in 
a spectrum which is highly symmetric about to = 
[lo|. This follows from the assumption that the vor- 
tex lattice causes the condensate to rotate as a rigid 
body. By contrast, the observed asymmetry in Fig. |21 
is due to the non-rigid rotation effects: If there are A^vtx 
vortices within the Thomas- Fermi radius ttf, the ve- 
locity of the condensate flow at ttf is approximately 
vtf ~ 2'!ThNvtx/MrT-p implying that the actual rotating 
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FIG. 3: Plots of angular momentum vs. energy uj„,m for 
(a) n = 0.5 (7 vortices) and (b) = 0.8 (31 vortices), and 
angular quantum number m vs. Un.m for (c) Q = 0.5 and (d) 
Q = 0.8. 



frequency of the condensate at the Thomas-Fermi radius 
is r^TF = WTF/27rrTF = fiN^tx/Mr^p. In the case of 7 
vortices, r^p = 5.9, which yields I^tf — 0.2. The rota- 
tion of the condensate on the edge of the cloud is there- 
fore much slower than the rotation of trap. One can show 
that the frequency of the circulating wave in the reference 
frame rotating at fixF is „ = ^^n,m + ™ 
and if we use this „ as the energy, then m vs. lo'^^ 
curve becomes almost symmetric about m = 0. With 
a sufficiently large number A'vtx (rigid body regime), 
ilxF ~ r2 T]. In the case of 31 vortices, the energy split- 
ting contribution proportional to to became rather small 
while the predicted 2^ splitting persisted. 

We note that for the positive to states the Lz vs. 0Jn,m 
curve is quite smooth, yet there are several exciations 
with lower angular momenta that appear as "dips" in the 
curve (see Fig. These dips occur at the same angular 
momentum even for different n. The quantum number 
TO for these states are approximately integral multiples of 
the number of outermost vortices, which is 6 for = 0.5 
and 12 for Vl ~ 0.8. This effect may be understood from 
the fact that the bulk excitations interact strongly with 
the vortex lattice when the intervortex spacing of the 
outermost vortices happens to be a multiple of the wave- 
length of the surface wave. We have observed that for 
these states, vortices, resonant with the surface wave, 
partly share the excitation energy with the surface wave 
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so as to reduce the angular momentum. This situation is 
reminiscent of the dispersion relation of electrons in an 
ion core lattice potential where, using a Kronnig-Penny 
model 01 one finds that deviations appears when an 
integral multiple electron wavefunction wavelengths co- 
incides with the inter-ion spacing. The analogy is not 
complete, however, because the vortices can be more eas- 
ily disturbed by the surface waves than can ion cores be 
displaced in a crystal. 

(a) (b) 



(c) (d) 




FIG. 4: Density plots of u for (a) the Tkachenko and (b) 
the hydrodynamic (n = 0, m = 1) mode; (c) and (d) are the 
corresponding plots of v. (e) gives the mean field potential 
that quasi particle u and v experience. 

Finally, to understand the Tkachenko and the hydro- 
dynamic shape modes at the level of quasiparticles, we 
plot the density of u and v for these two cases. Figures 
2Ia) and (b) show that, for the Tkachenko modes, u and 
V have localized peaks where the vortices are present, 
while for the hydrodynamic shape excitations, they are 
more evenly spread with depressions instead of peaks. 
The regions where the quasiparticles are localized are 
mainly excited. In particular, vortex precession is ob- 
served for the Tkachenko modes. Figure EJc) shows the 
effective mean field potential the quasiparticles experi- 
ence, Vtr + '2g'ip*'ip — fi, from the Bogoliubov equation, Eq. 
(O; the vortex cores create an array of effective poten- 
tial wells. It is clear that the Tkachenko modes may be 
viewed as the lower energy bounded states within those 
wells. The hydrodynamic shape modes may be viewed as 
"unbounded" states with higher energies - even though 
still bounded within the trapping potential. 

In summary, we have described new aspects of nor- 
mal modes of an atomic BEC supporting a vortex lat- 
tice. From a spectroscopic studies of the excitations, we 
found anomalous effects which happen due to the singu- 
larity of vorticity and due to the non-rigid rotation of the 
condensate. Some of our observations should be readily 
verifiable experimentally. Our results for the excitation 



spectrum is an ideal starting point for finite tempera- 
ture calculations related to the dissipation and decay of 
vortices. This will be the subject of future work. 

During the preparation of the manuscript we learned 
of a related preprint This work is supported by the 
NSF, ONR and ARO. L. O. Baks. is a Horton Fellow. 
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